The Physical Oceanography of Ice-Covered Moons.

In the outer solar system, a growing number of giant planet satellites are now known to be abodes for global oceans hidden below an outer layer of ice. These planetary oceans are a natural laboratory for studying physical oceanographic processes in settings that challenge traditional assumptions made for Earth's oceans. While some driving mechanisms are common to both systems, such as buoyancy-driven flows and tides, others, such as libration, precession, and electromagnetic pumping, are likely more significant for moons in orbit around a host planet. Here, we review these mechanisms and how they may operate across the solar system, including their implications for ice-ocean interactions. Future studies should continue to advance our understanding of each of these processes as well as how they may act together in concert. This interplay also has strong implications for habitability as well as testing oceanic hypotheses with future missions.

Melting of glacier ice enhanced by bursting air bubbles.

Feedbacks between ice melt, glacier flow and ocean circulation can rapidly accelerate ice loss at tidewater glaciers and alter projections of sea-level rise. At the core of these projections is a model for ice melt that neglects the fact that glacier ice contains pressurized bubbles of air due to its formation from compressed snow. Current model estimates can underpredict glacier melt at termini outside the region influenced by the subglacial discharge plume by a factor of 10-100 compared with observations. Here we use laboratory-scale experiments and theoretical arguments to show that the bursting of pressurized bubbles from glacier ice could be a source of this discrepancy. These bubbles eject air into the seawater, delivering additional buoyancy and impulses of turbulent kinetic energy to the boundary layer, accelerating ice melt. We show that real glacier ice melts 2.25 times faster than clear bubble-free ice when driven by natural convection in a laboratory setting. We extend these results to the geophysical scale to show how bubble dynamics contribute to ice melt from tidewater glaciers. Consequently, these results could increase the accuracy of modelled predictions of ice loss to better constrain sea-level rise projections globally.

PlasPi TDM: Augmentation of a low-cost camera platform for advanced underwater physical-ecological observations.

Marine ecosystem dynamics in the context of climate change is a growing scientific, political and social concern requiring regular monitoring through appropriate observational technologies and studies. Thus, a wide range of tools comprising chemical, biogeochemical, physical, and biological sensors, as well as other platforms exists for marine monitoring. However, their high acquisition and maintenance costs are often a major obstacle, especially in low-income developing countries. We designed an advanced low-cost synoptic marine ecosystem observation system that operates at relatively high temporal frequencies, named PlasPi TDM. This instrument is an improved version of the camera system (PlasPI marine cameras) developed in 2020 by Autun Purser from the Alfred Wegener Institute Helmholtz Center for Polar and Marine Research (Germany), and collaborators. It incorporates several innovative developments such as multispectral (records the spectrum of any object photographed), temperature and pressure sensors. The PlasPi TDM operates to a depth of 200 m. The various field deployments demonstrate the operational capability of the PlasPi TDM for different applications and illustrate its considerable potential for in-situ observations and marine surveillance in Africa. This device is intended as an open-source project and its continued development is encouraged for a more integrated, sustainable and low-cost ocean observing system.

Glacier retreat alters downstream fjord ecosystem structure and function in Greenland.

The melting of the Greenland Ice Sheet is accelerating, with glaciers shifting from marine to land termination and potential consequences for fjord ecosystems downstream. Monthly samples in 2016 in two fjords in southwest Greenland show that subglacial discharge from marine-terminating glaciers sustains high phytoplankton productivity that is dominated by diatoms and grazed by larger mesozooplankton throughout summer. In contrast, melting of land-terminating glaciers results in a fjord ecosystem dominated by bacteria, picophytoplankton and smaller zooplankton, which has only one-third of the annual productivity and half the CO2 uptake compared to the fjord downstream from marine-terminating glaciers.

Seasonality of the Meridional Overturning Circulation in the subpolar North Atlantic.

Understanding the variability of the Atlantic Meridional Overturning Circulation is essential for better predictions of our changing climate. Here we present an updated time series (August 2014 to June 2020) from the Overturning in the Subpolar North Atlantic Program. The 6-year time series allows us to observe the seasonality of the subpolar overturning and meridional heat and freshwater transports. The overturning peaks in late spring and reaches a minimum in early winter, with a peak-to-trough range of 9.0 Sv. The overturning seasonal timing can be explained by winter transformation and the export of dense water, modulated by a seasonally varying Ekman transport. Furthermore, over 55% of the total meridional freshwater transport variability can be explained by its seasonality, largely owing to overturning dynamics. Our results provide the first observational analysis of seasonality in the subpolar North Atlantic overturning and highlight its important contribution to the total overturning variability observed to date.

Crevasse refreezing and signatures of retreat observed at Kamb Ice Stream grounding zone.

Ice streams flowing into Ross Ice Shelf are presently responsible for around 10% of the mass flux from West Antarctica, with the noteworthy exception of Kamb Ice Stream, which stagnated in the late 1800s. The subsequent reduction in ice supply led to grounding-line retreat at the coastal margin where Kamb transitions into the floating Ross Ice Shelf. Grounding-line migration is linked to broader changes in ice-sheet mass balance and sea level, but our understanding of related ice, ocean and seafloor interactions is limited by the difficulty in accessing these remote regions. Here we report in situ observations from an underwater vehicle deployed at Kamb that show how fine-scale variability in ice and ocean structure combine to influence a diversity of ice-ocean interactions. We found a stratified water column within a tenth of a degree of freezing at the ice base and mapped basal crevasses with supercooled water and active marine ice formation. At the seafloor, we interpret parallel ridges as crevasse impressions left as the ice lifted off during grounding-line retreat. These observations from a recently ungrounded sub-shelf environment illuminate both the geomorphological signatures of past grounding-line retreat and the fine-scale sensitivity of ongoing ice-ocean interactions to ice topography.

Influence of Particle Size and Fragmentation on Large-Scale Microplastic Transport in the Mediterranean Sea.

Microplastic particles move three-dimensionally through the ocean, but modeling studies often do not consider size-dependent vertical transport processes. In addition, microplastic fragmentation in ocean environments remains poorly understood, despite fragments making up the majority of microplastic pollution in terms of the number of particles and despite its potential role in mass removal. Here, we first investigate the role of particle size and density on the large-scale transport of microplastics in the Mediterranean Sea and next analyze how fragmentation may affect transport and mass loss of plastics. For progressively smaller particle sizes, microplastics are shown to be less likely to be beached and more likely to reach open water. Smaller particles also generally get mixed deeper, resulting in lower near-surface concentrations of small particles despite their higher total abundance. Microplastic fragmentation is shown to be dominated by beach-based fragmentation, with ocean-based fragmentation processes likely having negligible influence. However, fragmentation remains a slow process acting on decadal time scales and as such likely does not have a major influence on the large-scale distribution of microplastics and mass loss over periods less than 3 years.

Oceanography dataset in Bonpies archipelago as remote island in Java Seas, Indonesia.

The oceanographic data presented in this article were collected in the Bonpies island, as part of the Java Seas. The data were collected by in-situ measurement from several expeditions from 2011 to 2016. The data presented here include bathymetry, sea surface temperature (SST), sea surface salinity (SSS), ocean transparency, dissolved oxygen, and pH. All data were collected by daily in-situ measurements in different seasons. There were 5768 data measured from several portable instruments, including echo-sounder and portable instruments. Calibration of the instrument was carried out before and after the survey in accordance with the protocol provided by the manufacturer. Measurements are made directly and recorded into a log sheet paper and the data format is Comma Separated Files (.csv). This in-situ data would be most useful for regional climate studies, including forecasting such as El Nino and Indian Ocean Dipole effects, oceanographic condition, and marine resources management.

Global Dynamics of the Stationary M2 Mode-1 Internal Tide.

A reduced-physics model is employed at 1/25° to 1/100° global resolution to determine (a) if linear dynamics can reproduce the observed low-mode M2 internal tide, (b) internal-tide sensitivity to bathymetry, stratification, surface tides, and dissipation parameterizations, and (c) the amount of power transferred to the nonstationary internal tide. The simulations predict 200 GW of mode-1 internal-tide generation, consistent with a general circulation model and semianalytical theory. Mode-1 energy is sensitive to damping, but a simulation using parameterizations for wave drag and wave-mean interaction predicts 84% of satellite observed sea-surface height amplitude variance on a 1° × 1° grid. The simulation energy balance indicates that 16% of stationary mode-1 energy is scattered to modes 2-4 and negligible energy propagates onto the shelves. The remaining 84% of energy is lost through parameterizations for high-mode scattering over rough topography (54%) and wave-mean interactions that transfer energy to the nonstationary internal tide (29%).

Right Place, Right Time: An Informal Memoir.

My career spanned the revolution in understanding of the large-scale fluid ocean, as modern electronics produced vast new capabilities. I started in the days of almost purely mechanical instruments operated by seagoing scientists, ones not so different from those used more than a century earlier. Elegant theories existed of hypothetical steady-state oceans. Today, we understand that the ocean is a highly turbulent fluid, interacting over global scales, and it is now studied by large teams using spacecraft and diverse sets of self-contained in situ instrumentation. Mine was an accidental career: I was lucky to be in the right place at the right time.

Significant contribution of small icebergs to the freshwater budget in Greenland fjords.

Icebergs represent nearly half of the mass loss from the Greenland Ice Sheet and provide a distributed source of freshwater along fjords which can alter fjord circulation, nutrient levels, and ultimately the Meridional Overturning Circulation. Here we present analyses of high resolution optical satellite imagery using convolutional neural networks to accurately delineate iceberg edges in two East Greenland fjords. We find that a significant portion of icebergs in fjords are comprised of small icebergs that were not detected in previously-available coarser resolution satellite images. We show that the preponderance of small icebergs results in high freshwater delivery, as well as a short life span of icebergs in fjords. We conclude that an inability to identify small icebergs leads to inaccurate frequency-size distribution of icebergs in Greenland fjords, an underestimation of iceberg area (specifically for small icebergs), and an overestimation of iceberg life span.

More meteorological events that drive compound coastal flooding are projected under climate change.

Compound flooding arises from storms causing concurrent extreme meteorological tides (that is the superposition of storm surge and waves) and precipitation. This flooding can severely affect densely populated low-lying coastal areas. Here, combining output from climate and ocean models, we analyse the concurrence probability of the meteorological conditions driving compound flooding. We show that, under a high emissions scenario, the concurrence probability would increase globally by more than 25% by 2100 compared to present. In latitudes above 40o north, compound flooding could become more than 2.5 times as frequent, in contrast to parts of the subtropics where it would weaken. Changes in extreme precipitation and meteorological tides account for most (77% and 20%, respectively) of the projected change in concurrence probability. The evolution of the dependence between precipitation and meteorological tide dominates the uncertainty in the projections. Our results indicate that not accounting for these effects in adaptation planning could leave coastal communities insufficiently protected against flooding.

VENM: An Algorithm to Accurately Calculate Neutral Slopes and Gradients.

Mesoscale eddies stir along the neutral plane, and the resulting neutral diffusion is a fundamental aspect of subgrid-scale tracer transport in ocean models. Calculating neutral diffusion traditionally involves calculating neutral slopes and three-dimensional tracer gradients. The calculation of the neutral slope traditionally occurs by computing the ratio of the horizontal to vertical locally referenced potential density derivative. However, this approach is problematic in regions of weak vertical stratification, prompting the use of a variety of ad hoc regularization methods that can lead to rather nonphysical dependencies for the resulting neutral tracer gradients. Here we use a VErtical Non-local Method "VENM," a search algorithm that requires no ad hoc regularization and significantly improves the numerical accuracy of calculating neutral slopes, neutral tracer gradients, and associated neutral diffusive fluxes. We compare and contrast VENM against a more traditional method, using an independent objective neutrality condition combined with estimates of spurious diffusion, heat transport, and water mass transformation rates. VENM is more accurate, both physically and numerically, and should form the basis for future efforts involving neutral diffusion calculations from observations and possibly numerical model simulations.

Unsupervised Learning Reveals Geography of Global Ocean Dynamical Regions.

Dynamically similar regions of the global ocean are identified using a barotropic vorticity (BV) framework from a 20-year mean of the Estimating the Circulation and Climate of the Ocean state estimate at 1° resolution. An unsupervised machine learning algorithm, K-means, objectively clusters the standardized BV equation, identifying five unambiguous regimes. Cluster 1 covers 43 ± 3.3% of the ocean area. Surface and bottom stress torque are balanced by the bottom pressure torque and the nonlinear torque. Cluster 2 covers 24.8 ± 1.2%, where the beta effect balances the bottom pressure torque. Cluster 3 covers 14.6 ± 1.0%, characterized by a "Quasi-Sverdrupian" regime where the beta effect is balanced by the wind and bottom stress term. The small region of Cluster 4 has baroclinic dynamics covering 6.9 ± 2.9% of the ocean. Cluster 5 occurs primarily in the Southern Ocean. Residual "dominantly nonlinear" regions highlight where the BV approach is inadequate, found in areas of rough topography in the Southern Ocean and along western boundaries.

River-discharge effects on United States Atlantic and Gulf coast sea-level changes.

Identifying physical processes responsible for historical coastal sea-level changes is important for anticipating future impacts. Recent studies sought to understand the drivers of interannual to multidecadal sea-level changes on the United States Atlantic and Gulf coasts. Ocean dynamics, terrestrial water storage, vertical land motion, and melting of land ice were highlighted as important mechanisms of sea-level change along this densely populated coast on these time scales. While known to exert an important control on coastal ocean circulation, variable river discharge has been absent from recent discussions of drivers of sea-level change. We update calculations from the 1970s, comparing annual river-discharge and coastal sea-level data along the Gulf of Maine, Mid-Atlantic Bight, South Atlantic Bight, and Gulf of Mexico during 1910-2017. We show that river-discharge and sea-level changes are significantly correlated ([Formula: see text]), such that sea level rises between 0.01 and 0.08 cm for a 1 [Formula: see text] annual river-discharge increase, depending on region. We formulate a theory that describes the relation between river-discharge and halosteric sea-level changes (i.e., changes in sea level related to salinity) as a function of river discharge, Earth's rotation, and density stratification. This theory correctly predicts the order of observed increment sea-level change per unit river-discharge anomaly, suggesting a causal relation. Our results have implications for remote sensing, climate modeling, interpreting Common Era proxy sea-level reconstructions, and projecting coastal flood risk.

Improved Visualization of Hydroacoustic Plumes Using the Split-Beam Aperture Coherence.

Natural seepage of methane into the oceans is considerable, and plays a role in the global carbon cycle. Estimating the amount of this greenhouse gas entering the water column is important in order to understand their environmental impact. In addition, leakage from man-made structures such as gas pipelines may have environmental and economical consequences and should be promptly detected. Split beam echo sounders (SBES) detect hydroacoustic plumes due to the significant contrast in acoustic impedance between water and free gas. SBES are also powerful tools for plume characterization, with the ability to provide absolute acoustic measurements, estimate bubble trajectories, and capture the frequency dependent response of bubbles. However, under challenging conditions such as deep water and considerable background noise, it can be difficult to detect the presence of gas seepage from the acoustic imagery alone. The spatial coherence of the wavefield measured across the split beam sectors, quantified by the coherence factor (CF), is a computationally simple, easily available quantity which complements the acoustic imagery and may ease the ability to automatically or visually detect bubbles in the water column. We demonstrate the benefits of CF processing using SBES data from the Hudson Canyon, acquired using the Simrad EK80 SBES. We observe that hydroacoustic plumes appear more clearly defined and are easier to detect in the CF imagery than in the acoustic backscatter images.

Application of the ideas and techniques of classical fluid mechanics to some problems in physical oceanography.

This review makes a case for describing many of the flows observed in our oceans, simply based on the Euler equation, with (piecewise) constant density and with suitable boundary conditions. The analyses start from the Euler and mass conservation equations, expressed in a rotating, spherical coordinate system (but the f-plane and ß-plane approximations are also mentioned); five examples are discussed. For three of them, a suitable non-dimensionalization is introduced, and a single small parameter is identified in each case. These three examples lead straightforwardly and directly to new results for: waves on the Pacific Equatorial Undercurrent (EUC) with a thermocline (in the f-plane); a nonlinear, three-dimensional model for EUC-type flows (in the ß-plane); and a detailed model for large gyres. The other two examples are exact solutions of the complete system: a flow which corresponds to the underlying structure of the Pacific EUC; and a flow based on the necessary requirement to use a non-conservative body force, which produces the type of flow observed in the Antarctic Circumpolar Current. (All these examples have been discussed in detail in the references cited.) This review concludes with a few comments on how these solutions can be extended and expanded.This article is part of the theme issue 'Nonlinear water waves'.

Locally stationary spatio-temporal interpolation of Argo profiling float data.

Argo floats measure seawater temperature and salinity in the upper 2000 m of the global ocean. Statistical analysis of the resulting spatio-temporal dataset is challenging owing to its non-stationary structure and large size. We propose mapping these data using locally stationary Gaussian process regression where covariance parameter estimation and spatio-temporal prediction are carried out in a moving-window fashion. This yields computationally tractable non-stationary anomaly fields without the need to explicitly model the non-stationary covariance structure. We also investigate Student t-distributed fine-scale variation as a means to account for non-Gaussian heavy tails in ocean temperature data. Cross-validation studies comparing the proposed approach with the existing state of the art demonstrate clear improvements in point predictions and show that accounting for the non-stationarity and non-Gaussianity is crucial for obtaining well-calibrated uncertainties. This approach also provides data-driven local estimates of the spatial and temporal dependence scales for the global ocean, which are of scientific interest in their own right.

Juvenile recruitment in loggerhead sea turtles linked to decadal changes in ocean circulation.

Given the threats of climate change, understanding the relationship of climate with long-term population dynamics is critical for wildlife conservation. Previous studies have linked decadal climate oscillations to indices of juvenile recruitment in loggerhead sea turtles (Caretta caretta), but without a clear understanding of mechanisms. Here, we explore the underlying processes that may explain these relationships. Using the eddy-resolving Ocean General Circulation Model for the Earth Simulator, we generate hatch-year trajectories for loggerhead turtles emanating from Japan over six decades (1950-2010). We find that the proximity of the high-velocity Kuroshio Current to the primary nesting areas in southern Japan is remarkably stable and that hatchling dispersal to oceanic habitats itself does not vary on decadal timescales. However, we observe a shift in latitudes of trajectories, consistent with the Pacific Decadal Oscillation (PDO). In a negative PDO phase, the Kuroshio Extension Current (KEC) is strong and acts as a physical barrier to the northward transport of neonates. As a result, hatch-year trajectories remain mostly below 35°N in the warm, unproductive region south of the Transition Zone Chlorophyll Front (TZCF). During a positive PDO phase, however, the KEC weakens facilitating the neonates to swim north of the TZCF into cooler and more productive waters. As a result, annual cohorts from negative PDO years may face a lack of resources, whereas cohorts from positive PDO years may find sufficient resources during their pivotal first year. These model outputs indicate that the ocean circulation dynamics, combined with navigational swimming behavior, may be a key factor in the observed decadal variability of sea turtle populations.

Cross-Shelf Exchange.

Cross-shelf exchange dominates the pathways and rates by which nutrients, biota, and materials on the continental shelf are delivered and removed. This follows because cross-shelf gradients of most properties are usually far greater than those in the alongshore direction. The resulting transports are limited by Earth's rotation, which inhibits flow from crossing isobaths. Thus, cross-shelf flows are generally weak compared with alongshore flows, and this leads to interesting observational issues. Cross-shelf flows are enabled by turbulent mixing processes, nonlinear processes (such as momentum advection), and time dependence. Thus, there is a wide range of possible effects that can allow these critical transports, and different natural settings are often governed by different combinations of processes. This review discusses examples of representative transport mechanisms and explores possible observational and theoretical paths to future progress.